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The present article reviews the recent developments in the physics of quantum turbulence. 
Quantum turbulence (QT) was discovered in superfluid 4 He in the 1950s, and the research 
has tended toward a new direction since the mid 90s. The similarities and differences between 
quantum and classical turbulence have become an important area of research. QT is comprised 
of quantized vortices that are definite topological defects, being expected to yield a model of 
turbulence that is much simpler than the classical model. The general introduction of the issue 
and a brief review on classical turbulence are followed by a description of the dynamics of 
quantized vortices. Then, we discuss the energy spectrum of QT at very low temperatures. At 
low wavenumbers, the energy is transferred through the Richardson cascade of quantized vor- 
tices, and the spectrum obeys the Kolmogorov law, which is the most important statistical law 
in turbulence; this classical region shows the similarity to conventional turbulence. At higher 
wavenumbers, the energy is transferred by the Kelvin-wave cascade on each vortex. This quan- 
tum regime depends strongly on the nature of each quantized vortex. The possible dissipation 
mechanism is discussed. Finally, important new experimental studies, which include investi- 
gations into temperature-dependent transition to QT, dissipation at very low temperatures, 
QT created by vibrating structures, and visualization of QT, are reviewed. The present article 
concludes with a brief look at QT in atomic Bose-Einstein condensates. 
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1. Introduction 

For high-velocity or high Reynolds number flows, the 
flow is generally turbulent. Turbulence has been investi- 
gated not only in basic science research, such as physics 
and mathematics research, but also in applied sciences, 
such as fluid engineering and aeronautics. Although tur- 
bulence has been studied intensely in a number of fields, 
it is still not yet well understood. This is chiefly because 
turbulence is a complicated dynamical phenomenon with 
strong nonlinearity that is quite different from an equi- 
librium state. Vortices may be key in understanding tur- 
bulence. For example, Leonardo da Vinci observed the 
turbulent flow of water and drew many sketches show- 
ing that turbulence had a structure comprised of vortices 
of different sizes. However, vortices are not well-defined 
for a typical classical fluid, and the relationship between 
turbulence and vortices remains unclear. 

Independently of these studies in classical fluid dynam- 
ics, turbulence in superfluid helium has been studied in 
the field of low-temperature physics. Liquid 4 He enters a 
superfluid state below the A point (2.17 K) with Bose- 
Einstein condensation of the 4 He atoms. ^ Characteristic 
phenomena of superfluidity were discovered experimen- 
tally in the 1930s by Kapitza et al. The hydrodynamics 
of superfluid helium is well described by the two-fluid 
model, for which the system consists of an inviscid su- 
perfluid (density p s ) and a viscous normal fluid (density 
p n ) with two independent velocity fields v s and v n . The 
mixing ratio of the two fluids depends on temperature. 
As the temperature is reduced below the A point, the ra- 
tio of the superfluid component increases, and the entire 



* E-mail: tsubota@sci.osaka-cu.ac.jp 

t Present address: Department of Physics, Osaka City University, 
Osaka 558-8585, Japan. 



fluid becomes a superfluid below approximately 1 K. The 
Bose-condensed system exhibits the macroscopic wave 
function ^(x, t) = |^(x, i)|e ie< - x ' t - ) as an order parameter. 
The superfluid velocity field is given by v s = (h/m)V9 
with boson mass m, representing the potential flow. Since 
the macroscopic wave function should be single-valued 
for the space coordinate x, the circulation F = § v ■ <M. 
for an arbitrary closed loop in the fluid is quantized by 
the quantum k = h/m. A vortex with such quantized 
circulation is called a quantized vortex. Any rotational 
motion of a superfluid is sustained only by quantized 
vortices. 

A quantized vortex is a topological defect characteris- 
tic of a Bose-Einstein condensate and is different from a 
vortex in a classical viscous fluid. First, the circulation 
is quantized, which is contrary to a classical vortex that 
can have any value of circulation. Second, a quantized 
vortex is a vortex of inviscid superflow. Thus, it cannot 
decay by the viscous diffusion of vorticity that occurs in 
a classical fluid. Third, the core of a quantized vortex is 
very thin, on the order of the coherence length, which is 
only a few angstroms in superfluid 4 He. Since the vortex 
core is very thin and does not decay by diffusion, it is 
always possible to identify the position of a quantized 
vortex in the fluid. These properties make a quantized 
vortex more stable and definite than a classical vortex. 

Early experimental studies on superfluid turbulence 
focused primarily on thermal counterflow, in which the 
normal fluid and superfluid flow in opposite directions. 
The flow is driven by an injected heat current, and it was 
found that the superflow becomes dissipative when the 
relative velocity between the two fluids exceeds a critical 
value. 2 ) Feynman proposed that this is a superfluid tur- 
bulent state consisting of a tangle of quantized vortices. 3 ) 
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Vinen later confirmed Feynman's findings experimentally 
by showing that the dissipation comes from the mutual 
friction between vortices and the normal flow. 4-6 ) Subse- 
quently several experimental studies have examined su- 
perfluid turbulence (ST) in thermal counterflow systems 
and have revealed a variety of physical phenomenon. 7 ) 
Since the dynamics of quantized vortices is nonlinear 
and non-local, it has not been easy to understand vor- 
tex dynamics observations quantitatively. Schwarz clar- 
ified the picture of ST consisting of tangled vortices by 
a numerical simulation of the quantized vortex filament 
model in the thermal counterflow. 8, 9 ) However, since the 
thermal counterflow has no analogy to conventional fluid 
dynamics, this study was not helpful in clarifying the 
relationship between ST and classical turbulence (CT). 
Superfluid turbulence is often called quantum turbulence 
(QT), which emphasizes the fact that it is comprised of 
quantized vortices. 

Comparing QT and CT reveals definite differences, 
which demonstrates the importance of studying QT. Tur- 
bulence in a classical viscous fluid appears to be com- 
prised of vortices, as pointed out by Da Vinci. However, 
these vortices are unstable, and appear and disappear re- 
peatedly. Moreover, the circulation is not conserved and 
is not identical for each vortex. Quantum turbulence con- 
sists of a tangle of quantized vortices that have the same 
conserved circulation. Thus, QT is an easier system to 
study than CT and has a much simpler model of turbu- 
lence than CT. 

Based on this consideration, QT research has tended 
toward a new direction since the mid 90s. A chief inter- 
est is to understand the relationship between QT and 
CT. 10 ' 11 ) The energy spectrum of fully developed clas- 
sical turbulence is known to obey the Kolmogorov law 
in an inertial range. The energy transfer in an inertial 
range is believed to be sustained by the Richardson cas- 
cade process, in which large eddies are broken up self- 
similarly to smaller eddies. Recent experimental and nu- 
merical studies support the Kolmogorov spectrum in QT. 
Another important problem is the dissipative mechanism 
in QT. At a finite temperature the mutual friction works 
as the dissipative mechanism. However, it is not so easy 
to understand what kind of mechanism causes dissipa- 
tion at very low temperatures in which the normal fluid 
component is negligible. 

This article reviews recent research on QT. After re- 
viewing classical turbulence in Section 2, Section 3 de- 
scribes the dynamics of quantized vortices, which is key 
in understanding QT. Section 4 summarizes the moti- 
vation and the problems of recent research on quantum 
turbulence. In Section 5, we discuss the issue of the en- 
ergy spectra and a possible dissipation mechanism at 
very low temperatures. Section 6 describes recent im- 
portant experimental studies, including research into the 
temperature-dependent transition to QT, dissipation at 
very low temperatures, QT created by vibrating struc- 
tures, and visualization of QT. Section 7 presents a sum- 
mary and briefly addresses QT in atomic Bose-Einstcin 
condensates. 
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2. Classical turbulence 

Before considering QT, we briefly review classical fluid 
dynamics and the statistical properties of CT. 12 ) Clas- 
sical viscous fluid dynamics is described by the Navier- 
Stokes equation: 

jU(x,i)+v(x,i)-Vv(x,i) = -ivP(x,t) + t/V 2 v(x,i), 

(1) 

where v(x, t) is the velocity of the fluid, P(x, t) is the 
pressure, p is the density of the fluid, and v is the kine- 
matic viscosity. The flow of this fluid can be character- 
ized by the ratio of the second term of the left-hand side 
of Eq. (1), hereinafter referred to as the inertial term, 
to the second term of the right-hand side, hereinafter 
called the viscous term. This ratio is the Reynolds num- 
ber R — vD/v, where v and D are the characteristic 
velocity of the flow and the characteristic scale, respec- 
tively. When v increases to allow the Reynolds number to 
exceed a critical value, the system changes from a lami- 
nar state to a turbulent state, in which the flow is highly 
complicated and contains many eddies. 

Such turbulent flow is known to show characteristic 
statistical behavior. 13 ' 14 ) We assume a steady state of 
fully developed turbulence of an incompressible classical 
fluid. The energy is injected into the fluid at a rate of e, 
the scale of which is comparable to the system size D in 
the energy-containing range. In the inertial range, this 
energy is transferred to smaller scales without being dis- 
sipated. In this range, the system is locally homogeneous 
and isotropic, which leads to the statistics of the energy 
spectrum known as the Kolmogorov law: 

E(k) = C£ 2 / 3 fc- 5 / 3 . (2) 

Here, the energy spectrum E(k) is defined as E = 
j dkE(k), where E is the kinetic energy per unit mass 
and k is the wavenumber from the Fourier transforma- 
tion of the velocity field. The spectrum of Eq. (2) is eas- 
ily derived by assuming that E(k) is locally determined 
by only the energy flux e and k. The energy transferred 
to smaller scales in the energy-dissipative range is dis- 
sipated at the Kolmogorov wavenumber kx = (e/V 3 ) 1 / 4 
through the viscosity of the fluid with dissipation rate e 
in Eq. (2), which is equal to the energy flux n in the iner- 
tial range. The Kolmogorov constant C is a dimensionless 
parameter of order unity. The Kolmogorov spectrum is 
confirmed experimentally and numerically in turbulence 
at high Reynolds numbers. The inertial range is thought 
to be sustained by the self-similar Richardson cascade 
in which large eddies are broken up into smaller eddies 
through many vortex reconnections. In CT, however, the 
Richardson cascade is not completely understood be- 
cause it is impossible to definitely identify each eddy. 
The Kolmogorov spectrum is based on the assumption 
that the turbulence is homogeneous and isotropic. How- 
ever, actual turbulence is not necessarily homogeneous 
or isotropic, and so the energy spectrum deviates from 
the Kolmogorov form. This phenomenon is called inter- 
mittency, and is an important problem in modern fluid 
dynamics. 12 -* Intermittcncy is closely related to the co- 
herent structure, which may be represented by vortices. 
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Research into QT could also shed light on this issue. 

3. Dynamics of quantized vortices 

As described in Section 1, most experimental studies 
on superfluid turbulence have examined thermal coun- 
terflow. However, the nonlinear and non-local dynamics 
of vortices have delayed progress in the microscopic un- 
derstanding of the vortex tangle. Schwarz overcame these 
difficulties 8 ' 9 ^ by developing a direct numerical simula- 
tion of vortex dynamics connected with dynamical scal- 
ing analysis, enabling the calculation of physical quanti- 
ties such as the vortex line density, anisotropic parame- 
ters, and mutual friction force. The observable quantities 
obtained by Schwarz agree well with some typical exper- 
imental results. 

Two formulations are generally available for studying 
the dynamics of quantized vortices. One is the vortex fil- 
ament model and the other is the Gross-Pitaevskii (GP) 
model. We will briefly review these two formulations. 

3. 1 Vortex filament model 

As described in Section 1 , a quantized vortex has quan- 
tized circulation. The vortex core is extremely thin, usu- 
ally much smaller than other characteristic scales in vor- 
tex motion. These properties allow a quantized vortex 
to be represented as a vortex filament. In classical fluid 
dynamics, 15 -* the vortex filament model is a convenient 
idealization. However, the vortex filament model is ac- 
curate and realistic for a quantized vortex in superfluid 
helium. 

The vortex filament formulation represents a quantized 
vortex as a filament passing through the fluid, having a 
definite direction corresponding to its vorticity. Except 
for the thin core region, the superflow velocity field has 
a classically well-defined meaning and can be described 
by ideal fluid dynamics. The velocity at a point r due to 
a filament is given by the Biot-Savart expression: 

(si — r) x efei 
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where k is the quantum of circulation. The filament is 
represented by the parametric form s = s(£, t) with the 
one-dimensional coordinate £ along the filament. The 
vector Si refers to a point on the filament, and the inte- 
gration is taken along the filament. Helmholtz's theorem 
for a perfect fluid states that the vortex moves with the 
superfluid velocity. Calculating the velocity w s at a point 
r = s on the filament causes the integral to diverge as 
Si — > s. To avoid this divergence, we separate the velocity 
s of the filament at the point s into two components: 8 - 1 
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The first term is the localized induction field arising from 
a curved line element acting on itself, and £ + and i- are 
the lengths of the two adjacent line elements after dis- 
cretization, separated by the point s. The prime denotes 
differentiation with respect to the arc length £. The mu- 
tually perpendicular vectors s', s", and s' x s" point along 
the tangent, the principal normal, and the binormal, re- 



spectively, at the point s, and their respective magni- 
tudes are 1, i£ _1 , and R~ x , with the local radius R of 
curvature. The parameter ao is the cutoff correspond- 
ing to the core radius. Thus, the first term represents 
the tendency to move the local point s along the binor- 
mal direction with a velocity inversely proportional to R. 
The second term represents the non-local field obtained 
by integrating the integral of Eq. (3) along the rest of 
the filament, except in the neighborhood of s. 

Neglecting the non-local terms and replacing Eq. (4) 
by s = /3s' x s" is referred to as the localized induc- 
tion approximation (LIA). Here, the coefficient j3 is de- 
fined by (3 = (k/47t) In (c(R)/a ), where c is a constant 
of order 1 and (l+l-) 1 / 2 is replaced by the mean ra- 
dius of curvature (R) along the length of the filament. 
This approximation is believed to be effective for analyz- 
ing isotropic dense tangles due to cancellations between 
non-local contributions. 

A better understanding of vortices in a real system 
results is obtained when boundaries are included in the 
analysis. For this purpose, a boundary-induced velocity 
field v Sj b is added to v s , so that the superflow can satisfy 
the boundary condition of an inviscid flow, that is, that 
the normal component of the velocity should disappear 
at the boundaries. To allow for another, presently un- 
specified, applied field, we include v Sjfl . Hence, the total 
velocity s of the vortex filament without dissipation is 
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At finite temperatures, it is necessary to take into ac- 
count the mutual friction between the vortex core and 
the normal flow v„. Including this term, the velocity of 
s is given by 

s = s + as' x (v„ - s ) - a's' x [s' x (v„ - s )], (6) 

where a and a' are temperature-dependent friction coef- 
ficients, 8 ) and so is calculated from Eq. (5). 

The numerical simulation method based on this model 
is described in detail elsewhere. 8,9 ' A vortex filament 
is represented by a single string of points separated by 
a distance A£. The vortex configuration at a given time 
determines the velocity field in the fluid, thus moving the 
vortex filaments according to Eqs. (5) and (6). Vortex re- 
connection should be properly included when simulating 
vortex dynamics. A numerical study of a classical fluid 
shows that the close interaction of two vortices leads to 
their reconnection, chiefly because of the viscous diffu- 
sion of the vorticity. 17 ) Schwarz assumed that two vor- 
tex filaments reconnect when they come within a critical 
distance of one another and showed that statistical quan- 
tities such as vortex line density were not sensitive to how 
the reconnections occur. 8,9 ) Even after Schwarz's study, 
it remained unclear as to whether quantized vortices can 
actually reconnect. However, Koplik and Levine solved 
directly the GP equation to show that two closely quan- 
tized vortices reconnect even in an inviscid superfluid. 18 ) 
More recent simulations have shown that reconnections 
are accompanied by emissions of sound waves having 
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wavelengths on the order of the healing length. 19 ' 20 ) 

Starting with several remnant vortices under thermal 
counterflow, Schwarz studied numerically how these vor- 
tices developed into a statistical steady vortex tangle. 9 ) 
The tangle was self-sustained by the competition be- 
tween the excitation due to the applied flow and the 
dissipation through the mutual friction. The numerical 
results were quantitatively consistent with typical experi- 
mental results. This was a significant accomplishment for 
numerical research. 

Here, we shall introduce some quantities that are char- 
acteristic of a vortex tangle. The line length density L 
is defined as the total length of vortex cores in a unit 
volume. The mean spacing I between vortices is given by 
i = L~ 1 / 2 . 

3.2 The Gross- Pitaevskii (GP) model 

In a weakly interacting Bose system, the macroscopic 
wave function ^(x, t) appears as the order parame- 
ter of Bose-Einstein condensation, obeying the Gross- 
Pitaevskii (GP) equation: 21 ) 

Here, g — Ai:h 2 m/a represents the strength of the inter- 
action characterized by the s-wave scattering length a, m 
is the mass of each particle, and n is the chemical poten- 
tial. Writing "J = \ty\ cxp(i8), the squared amplitude |W| 2 
is the condensate density and the gradient of the phase 
8 gives the superfluid velocity v s = (h/m)V9, which is a 
frictionless flow of the condensate. This relation causes 
quantized vortices to appear with quantized circulation. 
The only characteristic scale of the GP model is the co- 
herence length defined by £ = hj ( v / 2mg|'i , |), which gives 
the vortex core size. The GP model can explain not only 
the vortex dynamics but also phenomena related to vor- 
tex cores, such as reconnection and nucleation. However, 
the GP equation is not applicable quantitatively to su- 
perfluid 4 He, which is not a weakly interacting Bose sys- 
tem. It is, however, applicable to Bose-Einstein conden- 
sation of a dilute atomic Bose gas. 21 ) 

4. Modern research trends in QT 

Most older experimental studies on QT were devoted 
to thermal counterflow. Since this flow has no classical 
analogue, these studies do not contribute greatly to the 
understanding of the relationship between CT and QT. 
Since the mid 90s, important experimental studies on 
QT that did not focus on thermal counterflow have been 
published. 

4-.1 New experiments on energy spectra 

The first important contribution was made by Maurer 
and Tabeling, 22 ) who confirmed experimentally the Kol- 
mogorov spectrum in superfluid 4 He for the first time. A 
turbulent flow was produced in a cylinder by driving two 
counter-rotating disks. The authors observed the local 
pressure fluctuations to obtain the energy spectrum. The 
experiments were conducted at three different tempera- 
tures 2.3 K, 2.08 K, and 1.4 K. Both above and below the 
A point, the Kolmogorov spectrum was confirmed. The 
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observed behavior above the A point is not surprising be- 
cause the system is a classical viscous fluid. However, it 
is not trivial to understand the Kolmogorov spectrum at 
two different temperatures below the A point. 

The next significant step was a series of experiments 
on grid turbulence performed for superfluid 4 He above 1 
K by the Oregon group. 23 ~ 27 ) The flow through a grid 
is usually used to generate turbulence in classical fluid 
dynamics. 12 ' At a sufficient distance behind the grid, 
the flow displays a form of homogeneous isotropic turbu- 
lence. This method has also been applied to superfluid 
helium. In the experiments by the Oregon group, the 
helium was contained in a channel with a square cross 
section, through which a grid was pulled at a constant 
velocity. A pair of second-sound transducers was set into 
the walls of the channel. When a vortex tangle appeared 
in a channel, it was detected by second-sound attenua- 
tion. 10 ) The decay of the vorticity of the tangle created 
behind the towed grid was observed by the pair of trans- 
ducers. In combining the observations with the decay of 
the turbulence, the authors made some assumptions. In 
fully developed turbulence, the energy dissipation rate 
can be shown to be given as e = f(w 2 ), where (lj 2 ) is 
the mean square vorticity (rot v) in the flow. 28 ) The au- 
thors assumed that a similar formula applies in superfluid 
helium above 1 K. They noted that the quantity k 2 L 2 
would be a measure of the mean square vorticity in the 
superfluid component. Hence, they assumed that in grid 
turbulence the dissipation rate is given by e = v'k 2 L 2 
with an effective kinematic viscosity v' . In order to com- 
bine this representation with the observations of second- 
sound attenuation for grid turbulence, the authors fur- 
thermore assumed that a quasi-classical flow appears at 
length scales much greater than I. The flow is thought 
to come from a mechanism coupling the superfluid and 
the normal fluid by mutual friction, causing the fluid to 
behave like a one-component fluid. 29 ) By choosing suit- 
able values of v' as a function of temperature, 27 ) it was 
found that kL decays as t~ 3 / 2 . 

This characteristic decay L oc t~ 3 / 2 is quite impor- 
tant, because it is related to the Kolmogorov spectrum. 
Thus, we herein present the simple argument given in a 
previous review article. 10 ) We assume that a turbulent 
fluid obeys the Kolmogorov spectrum of Eq. (2) in the 
inertial range of D~ x < k < kd with D^ 1 « kd- The 
total energy is approximately given by 

E = Ce 2 ' 3 ^ 5 ' 3 « ^Ce 2 ' 3 D 2 ' 3 . (8) 
Jd- 1 2 

If the turbulence decays slowly with time and the dissi- 
pation rate e is assumed to be time-dependent, we can 
write 

e=-^L = -Ce- 1 / 3 D 2 / 3 -. (9) 
dt dt w 

The solution gives the time-dependence of e: 

e = 27C 3 D 2 (t + t )- 3 , (10) 
where t is a constant. Combining Eq. (10) with the 
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above formula of e gives the decay of L: 

This behavior has been observed, and a quantitative 
comparison with observations at any temperature gives 
v' as a function of temperature. The observation of the 
decay L ~ i~ 3 / 2 indicates that the Kolmogorov spectrum 
applies in turbulence, although it is not necessarily direct 
proof. Note that this simple analysis is applicable only 
when the maximum length scale of the turbulent energy 
saturates at the size of the channel. For the complete dy- 
namics, a more complicated decay of vortices has been 
observed and has been found to be consistent with the 
classical model of the Kolmogorov spectrum. 24 ~ 26 ' This 
type of decay is also observed at very low temperatures 
in the turbulence induced by an impulsive spin down for 
supcrffuid 4 He 30 ' and in grid turbulence for supcrfiuid 
3 He-B. 31 ' 32 ' 

4-2 Energy spectra at finite temperatures 

These observations lead us to inquire as to the na- 
ture of the velocity field that gives rise to the observed 
energy spectrum. Vincn considered the situation theoret- 
ically 29 ' and proposed that it is likely that the superfluid 
and normal fluid are coupled by mutual friction at scales 
larger than the characteristic scale I of the vortex tangle. 
If so, the two fluids behave as a one-component fluid at 
these scales, where mutual friction does not cause dis- 
sipation. Since the normal fluid is viscous, the two cou- 
pled fluids can be turbulent and obey the Kolmogorov 
spectrum. The observation of the energy spectrum by 
Maurer et al. and Stalp et al. should support the idea 
of the coupled dynamics of two fluids. This behavior has 
been confirmed numerically as well. 33 ' At small scales, 
the two fluids should be decoupled, so that both mutual 
friction and normal fluid viscosity operate. What then 
happens to the energy spectrum? Although theoretical 
consideration has been given to this problem, 34 ' the an- 
swer remains controversial and has not yet been clarified. 
While this is an important problem, it is not investigated 
in the present study. 

4-3 Modern research trends 

The following three trends are currently the primary 
research areas in QT. The first is the energy spectra 
and the dissipation mechanism at zero temperature. 35 ' 
The second is QT created by vibrating structures. 36 ' The 
third is visualization of QT. 37 ' The remainder of this ar- 
ticle is devoted to the review of these topics. 

5. Energy spectra and dissipation at zero tem- 
perature 

What happens to QT at zero temperature is not so 
trivial. The first problem is determining the nature of 
the energy spectrum of turbulence for the pure supcrfiuid 
component. 35 ' If QT has a classical analogue, the energy 
spectrum is expected to obey the Kolmogorov law and 
have an inertial range in which the energy is transferred 
self-similarly from large to small scales. In QT at zero 
temperature, any rotational motion should be carried 



by quantized vortices. Since quantized vortices are defi- 
nite topological defects, the cascade can be attributed di- 
rectly to their dynamics, which is different from the case 
for CT. The second problem is determining how the en- 
ergy is transferred from large to small scales. There is no 
dissipative mechanism at large scales. However, some dis- 
sipative mechanism should operate at small scales, as de- 
scribed below. A vortex tangle has a characteristic scale 
£, which is defined by the mean spacing between vor- 
tex lines. At scales greater than £, a Richardson cascade 
should transfer the energy through the breakup of vor- 
tices. However, the Richardson cascade becomes ineffec- 
tive at small scales, especially below t. What mechanism 
cascades the energy instead of the Richardson cascade 
at these scales? The third problem is understanding the 
dissipation in the system. The first possibility is acoustic 
emission at vortex reconnections. In classical fluid dy- 
namics, vortex reconnections cause acoustic emission. In 
quantum fluids, numerical simulations of the GP model 
shows acoustic emission at every reconnection event. 19 ' 
However, this mechanism is thought to be unimportant 
because of the very short coherence length. The second 
possible mechanism is the radiation of sound (phonons) 
by the oscillatory motion of vortex cores. We will return 
to these mechanisms, as well as other possibilities. 

5.1 Energy spectra 

No experimental studies have addressed this issue di- 
rectly, although a few numerical studies have been con- 
ducted. The first study was performed by Nore et al. 
using the GP model. 38,39 ' They solved the GP equation 
numerically starting from Taylor-Green vortices, and fol- 
lowed the time development. The quantized vortices be- 
come tangled and the energy spectra of the incompress- 
ible kinetic energy seemed to obey the Kolmogorov law 
for a short period, although the energy spectra even- 
tually deviated from the Kolmogorov law. The second 
study was performed by the vortex filament model, 40 ' 
and the third study was performed by the modified GP 
model. 41 ' 42 ' 

5.1.1 Energy spectra by the vortex filament model 

Using the vortex filament model, Araki et al. generated 
a vortex tangle arising from Taylor-Green vortices and 
obtained an energy spectrum consistent with the Kol- 
mogorov law. 40 ' It would be informative to describe how 
to obtain the energy spectra under the vortex filament 
model. The energy spectrum is calculated by Fourier 
transform of the superfluid velocity v s (r) , which is deter- 
mined by the configuration of quantized vortices. Thus, 
the energy spectrum can be calculated directly from the 
configuration of the vortices. Using the Fourier transform 
v s (k) = (2n)~ 3 f drv s (r) exp(— ik-r) and Parseval's the- 
orem / dk|v s (k)| 2 = (2n)~ 3 f dr\v s (r)\ 2 , the kinetic en- 
ergy of the superfluid velocity per unit mass is 

E=\J dr|v s (r)| 2 = J dk|v s (k)| 2 . (12) 

The vorticity ui(r) = rotv s (r) is represented in Fourier 
space as v s (k) = ik x u(k)/|k| 2 , so that we have E — 
((2tt) 3 /2) / dk|w(k)| 2 /|k| 2 . The vorticity is concentrated 
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on the vortex filament, represented by 

w(r) = k [ dts'(06(s(0 - r), (13) 



which can be rewritten as 



W(k) = (2^ / d€s ' (?) exp Hs(0 ' k )- ( 14 ) 

Using the definition of the energy spectrum -B(fc) from 
E = J Q dkE(k), these relations yield 



^1^2S'(6) • S'(6) 



xex P Hk-(s(a)-s(6))), (15) 

where dil^ = k 2 s,\n6kd6]~d(l)k is the volume element in 
spherical coordinates. 

Starting from a Taylor-Green vortex and following the 
vortex motion under the full Biot-Savart law without 
mutual friction, Araki et al. obtained a roughly homoge- 
neous and isotropic vortex tangle. 40 ) This was a decay- 
ing turbulence, dissipated by the cutoff of the smallest 
vortices, the size of which is comparable to the space res- 
olution. Initially, the energy spectrum has a large peak 
at the largest scale, where the energy is concentrated. 
The spectrum changes as the vortices become homoge- 
neous and isotropic. The energy spectrum of the vortex 
tangle at some late stage was quantitatively consistent 
with the Kolmogorov spectrum in the small k region. By 
monitoring the development of the vortex size distribu- 
tion, the decay of a tangle is found to be sustained by 
a Richardson cascade process. These results support the 
quasi-classical picture of the inertial range in QT at very 
low temperatures. 

5.1.2 Energy spectra by the GP model 

As the third trial, the Kolmogorov spectra was con- 
firmed for both decaying 41 ) and steady 42 ) QT by the 
modified GP model. The normalized GP equation is 



[-V 2 - M + ff |0>(x,i)| 2 ]<i>(x,i), (16) 



which determines the dynamics of the macroscopic wave 
function <E>(x, t) = /(x, t) exp[i</>(x, t)\. The condensate 
density is |$(x, t)\ 2 — /(x, t) 2 , and the supcrfluid veloc- 
ity v(x, t) is given by v(x, t) = 2V</>(x, t). The vortic- 
ity w(x, t) = rotv(x, t) vanishes everywhere in a single- 
connected region of the fluid and thus all rotational flow 
is carried by quantized vortices. In the core of each vor- 
tex, $(x, t) vanishes so that the circulation around the 
core is quantized by 47r. The vortex core size is given by 
the healing length £ = 1/ f^/g. 

Note that the hydrodynamics described by the GP 
model is compressible, which is different from the case 
of the vortex filament model. The total energy of the GP 
model 



E(t) 



dx$*(x, t) 



-v- 



: + |/(x,t) 2 



*(x,t), (17) 



is represented by the sum of the interaction energy 
Eint(t), the quantum energy E q (t), and the kinetic en- 



ergy E km (t): 38 ' 39) 

E in t(t) = ^ f d X f(x,t) 4 ,E q (t) = jj I r/x[V/(x./)p 



E kin {t) = ^J dx[/(x,i)V0(x,t)] 2 . (18) 



The kinetic energy is furthermore divided into a com- 
pressible part E kin (t) due to compressible excitations 
and an incompressible part E l kin (t) due to vortices. If the 
Kolmogorov spectrum is observed, the spectrum should 
be that for the incompressible kinetic energy E l kin {t). 

The failure to obtain the Kolmogorov law under the 
pure GP model 38 ' 39 ) would be attributable to the fol- 
lowing reasons. Note that the situation here is decay- 
ing turbulence. Although the total energy E{t) was con- 
served, E l kin (t) decreased with increasing E kin {t). This 
was because many compressible excitations were cre- 
ated through repeated vortex reconnections 19 ' 20 ) and 
disturbed the Richardson cascade of quantized vortices 
even at large scales. 

Kobayashi and Tsubota overcame the difficulties of 
Nore et al. and obtained the Kolmogorov spectra in QT 
and clearly revealed the energy cascade. 41,42 ) They per- 
formed numerical calculation for the Fourier transformed 
GP equation with dissipation: 



(*-7(k))^l>(k,t) = [fc 2 -^)]*(k,t) 
+7F2 Yl *)**(k 2 , t) x *(k - ki + k 2 , t). (19) 



ki,k 2 



Here, <E>(k, t) is the spatial Fourier component of $(x, t) 
and V is the system volume. The healing length is given 
by £ = 1/| &\^/g. The dissipation should have the form 
7(k) = 7o#(fc — 27r/£) with the step function 9, which 
dissipates only the excitations smaller than £. This form 
of dissipation can be justified by the coupled analysis of 
the GP equation for the macroscopic wave function and 
the Bogoliubov-de Gennes equations for thermal excita- 
tions. 43 ) 

First, Kobayashi et al. confirmed the Kolmogorov 
spectra for decaying turbulence. 41 ) To obtain a turbu- 
lent state, they started the calculation from an initial 
configuration in which the density was uniform and the 
phase of the wave function had a random spatial distribu- 
tion. The initial wave function was dynamically unstable 
and soon developed into fully developed turbulence with 
many quantized vortex loops. The spectrum E l kin {k,t) 
of the incompressible kinetic energy was then found to 
obey the Kolmogorov law. 

A more elaborate analysis of steady QT was performed 
by introducing energy injection at large scales as well 
as energy dissipation at small scales. 42 ) Energy injection 
at large scales was effected by moving a random poten- 
tial V^x, t). Numerically, Kobayashi et al. placed random 
numbers between and Vq in space-time (x, t) at inter- 
vals of Xq for space and To for time and connected them 
smoothly using a four-dimensional spline interpolation. 
The moving random potential exhibited a Gaussian two- 
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(x-x') 2 {t-t'f 



2X1 



2T§ 



point correlation: 
(V(x,t)V(x',t')) = ^ 2 exp 

(20) 

This moving random potential had a characteristic spa- 
tial scale of X . Small vortex loops were first nucleated 
by the random potential, growing to the scale of X 
by its motion subjected to Eq. (20). The vortex loops 
were then cast into the Richardson cascade. If steady 
QT is obtained by the balance between the energy in- 
jection and the dissipation, it should have an energy- 
containing range of k < 2-k/Xo, an inertial range of 
2tt/X < k < 27r/£, and an energy-dissipative range of 
2tt/£ < k. 




0* 





1 



2 3 
t 



25 26 27 28 29 30 
t 



(a) 



(b) 



Fig. 1. Time development of E(t), E kin (t), E| in (t), and E l kin (t) 
at (a) the initial stage < t < 5 and (b) a later stage 25 < t < 
30. 42 ) 



incompressible kinetic energy E l kin (t) due to vortices is 
dominant in the total kinetic energy Eki n (t). The four 
energies are almost constant for t > 25, and steady QT 
was obtained. 

Such a steady QT enables us to investigate the energy 
cascade. Here, we expect an energy flow in wavenum- 
ber space similar to that in Fig. 2. The upper half of 
the diagram shows the kinetic energy .EL (t) of quan- 
tized vortices, and the lower half shows the kinetic en- 
ergy E kin {t) of compressible excitations. In the energy- 
containing range k < 2it/Xq, the system receives in- 
compressible kinetic energy from the moving random po- 
tential. During the Richardson cascade process of quan- 
tized vortices, the energy flows from small to large k in 
the inertial range 27r/Xo < k < 27r/£. In the energy- 
dissipative range 27r/£ < k, the incompressible kinetic 
energy transforms to compressible kinetic energy through 
reconnections of vortices or the disappearance of small 
vortex loops. The moving random potential also creates 
long-wavelength compressible sound waves, which are an- 
other source of compressible kinetic energy and also pro- 
duce an interaction with vortices. However, the effect of 
sound waves is weak because E l kin (t) is much larger than 
E kin (t), as shown in Fig. 1. 

This cascade can be confirmed quantitatively by check- 
ing whether the energy dissipation rate e of E|, in (t) is 
comparable to the flux of energy H(k,t) through the 
Richardson cascade in the inertial range. Although the 
details are described in Reference 42, H(k, t) is found to 
be approximately independent of k and comparable to 
s. As shown in Fig. 3, the energy spectrum is quantita- 
tively consistent with the Kolmogorov law in the inertial 
range 2ir/X < k < 2tt/£_ (b), which is equivalent to 
0.79 < k < 6.3. 



27t/Xo 



energy-containing 
range 



inert ial-range 



271/5 

energy-dis s ipative 
range 




Decay of 
| vortices through 
reconnections J 



Fig. 2. Flow of the incompressible kinetic energy Ei- (t) (upper 
half of diagram) and compressible kinetic energy E?. (t) (lower 
half) in wavenumber space. 42 ) 
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Fig. 3. (a) A typical vortex tangle, (b) Energy spectrum 
Ei . (k,t) for QT. The plotted points are from an ensemble av- 
erage of 50 randomly selected states at t > 25. The solid line is 
the Kolmogorov law. 42 ) 



A typical simulation of steady turbulence was per- 
formed for V = 50, X a = 4, and T = 6.4 x 10~ 2 . The 
dynamics started from the uniform wave function. Fig- 
ure 1 shows the time development of each energy com- 
ponent. The moving random potential nucleates sound 
waves as well as vortices, but both figures show that the 



Kobayashi et al. studied the decay of QT under the 
same formulation. 44 ) After obtaining a steady tangle, 
they switched off the motion of the random potential 
and found that L decayed as £~ 3 / 2 . 
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5.2 The Kelvin-wave cascade 

The arguments in the last section were chiefly limited 
to the large scale, usually larger than the mean spacing 
£ of a vortex tangle, in which the Richardson cascade is 
effective for transferring the energy from large to small 
scales. Here, we should consider what happens at smaller 
scales, for which the Richardson cascade should be less 
effective. The most probable scenario is the Kelvin-wave 
cascade. A Kelvin- wave is a deformation of a vortex line 
into a helix with the deformation propagating as a wave 
along the vortex line. 45 ) Kelvin-waves were first observed 
by making torsional oscillations in uniformly rotating su- 
perfluid 4 He. 46 ' 47 -' The approximate dispersion relation 
for a rectilinear vortex is ujk = («;fc 2 )/(47r)(ln(l/fcao) +c) 
with a constant c <~ 1. Note that this k is the wavenum- 
ber of an excited Kelvin-wave and is different from the 
wavenumber used for the energy spectrum in the last 
subsection. At a finite temperature, a significant fraction 
of normal fluid damps Kelvin- waves through mutual fric- 
tion. At very low temperatures, however, mutual friction 
does not occur, and the only possible mechanism of dis- 
sipation is the radiation of phonons. 48 ) Phonon radia- 
tion becomes effective only when the frequency becomes 
very high, typically on the order of GHz (k <~ 10 _1 
nm _1 ), so a mechanism is required to transfer the en- 
ergy to such high wavenumbers in order for Kelvin-waves 
to be damped. An early numerical simulation based on 
the vortex filament model showed that Kelvin- waves are 
unstable to the buildup of side bands. 49 ) This indicates 
the possibility that nonlinear interactions between differ- 
ent Kelvin- wave numbers can transfer energy from small 
to large wavenumbers, namely the Kelvin-wave cascade. 
This idea was first suggested by Svistunov 50 ) and was 
later developed and confirmed through theoretical and 
numerical analyses by Vinen, Tsubota and Mitani, 51 ) 
and Kozik and Svistunov. 52-54 ) 

Observation of the Kelvin-wave cascade is crucial. 
Such studies are not easy for a vortex tangle. The easiest 
approach would be to consider rotation. In a rotating ves- 
sel, quantized vortices form a vortex lattice parallel to the 
rotational axis. By oscillating the vortices, Kelvin-waves 
can be excited. This approach was used in the pioneer- 
ing experiments of Hall. 46,47 ' The challenge is detecting 
the Kelvin-wave cascade when it occurs. There are two 
possible methods. The first is the direct visualization of 
vortices. Recently, remarkable progress has been made in 
visualizing superflow and quantized vortices, 55 ' 56 ' as de- 
scribed in Section 6.4. Bewley et al. visualized quantized 
vortices by trapping micron-sized solid hydrogen parti- 
cles. 56 ) They also observed a vortex array under rota- 
tion. Applying this technique to the Kelvin-wave cascade 
could reveal important information directly. The second 
method is to observe acoustic emission resulting from 
a Kelvin-wave cascade. Since the frequency of emitted 
phonons is estimated to be of GHz order, this observation 
may not be easy and presents a challenging experimental 
problem. 

5.3 Classical- quantum crossover 

An important trend is investigation of the nature of 
the transition between Richardson (classical) and Kelvin- 



wave (quantum) cascades. Several theoretical considera- 
tions on this topic have been investigated, 57 ' 58 ^ although 
few numerical or experimental studies have been con- 
ducted. 

L'vov et al. theoretically discussed a bottleneck 
crossover between the two regions. 57 ) Their investiga- 
tion was based on the idea that a vortex tangle obeying 
the Kolmogorov spectrum is polarized to some degree, 
and they viewed the vortex tangle as a set of vortex 
bundles. When they tried to connect the energy spec- 
tra between two cascades, a serious mismatch occurred 
at the crossover scale between the classical and quantum 
spectra. This mismatch prevented the energy flux from 
propagating fully through the crossover region, which 
is referred to as the bottleneck effect. In order to rem- 
edy this mismatch, L'vov et al. used warm cascade solu- 
tions 59 ) and proposed a thermal-equilibrium type spec- 
trum between the classical and quantum spectra. The 
analysis of L'vov et al. was based on the assumption that 
the coarse-grained macroscopic picture of quantized vor- 
tices remains valid down to a scale of £. Without this 
assumption, Kozik and Svistunov theoretically investi- 
gated the details of the structure of the vortex bundle in 
the crossover region. 58 ) Depending on the types of vor- 
tex reconnections, the crossover range near £ was further 
divided into three subranges, resolving the mismatch be- 
tween two cascades. However, this topic is not yet fixed, 
and remains controversial. 

5.4 Possible dissipation mechanism at zero temperature 
Here, we should investigate the possible dissipation 
mechanisms of QT at very low temperatures, in which 
the normal fluid component is so negligible that mutual 
friction does not occur. In this case, there is no dissipa- 
tion at relatively large scales, and dissipation can only 
occur at small scales. Energy at large scales should be 
transferred to smaller scales by a Richardson cascade 
and then a Kelvin- wave cascade until the dissipation be- 
comes effective. The description presented herein is only 
schematic. Please refer to other papers for a more de- 
tailed discussion. 10 ' 29 ) 

Dissipation of QT at zero temperature was first pro- 
posed by Feynman. 3 ) Feynman suggested that a large 
vortex loop should be broken up into smaller loops 
through repeated vortex reconnections, which is essen- 
tially identical to the Richardson cascade. Feynman 
thought that the smallest vortex ring having a radius 
comparable to the atomic scale must be a roton, al- 
though this belief is not currently accepted universally. 
Later, Vinen considered the decay of superfluid turbu- 
lence in order to understand his experimental results for 
thermal counterflow, 4 ' 5 ) leading him to propose Vinen's 
equation. 6 ^ Although the experiments were performed at 
finite temperatures, the physics is important in order 
to study the decay of QT, as described briefly herein. 
Vortices in turbulence are assumed to be approximately 
evenly spaced with a mean spacing of £ = L~ x / 2 . The 
energy of the vortices then spreads from the vortices of 
wavenumber l/£ to a wide range of wavenumbers. The 
overall decay of the energy is governed by the character- 
istic velocity v s = k/2it£ and the time constant £/v s of 
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the vortices of size £, giving 



dt 



-X. 



i/v s 



-x l , 



(21) 



where \ is a dimensionless parameter that is generally 
dependent on temperature. Rewriting this equation using 
L, we obtain 



dL k 2 

~dt = ~ X 2^ ' 



(22) 



This is called Vinen's equation. The original Vinen's 
equation also includes a term of vortex amplification, 
which is not shown here. This equation describes the de- 
cay of L, the solution of which is 



1 1 K 



(23) 



where Lq is L at t — 0. The thermal counterflow obser- 
vations are well described by this solution, which gives 
the values of x as a function of temperature. 6 ) The de- 
cay at finite temperatures is due to mutual friction. The 
value of x at zero temperature is obtained by a numerical 
simulation of the vortex filament model. 16 ) 

What causes the different types of decay of L ~ i~ 3 / 2 
of Eq. (11) and L ~ i" 1 of Eq. (23)? The different types 
of decay originate from the different structures of vor- 
tex tangles in QT. When the energy spectrum of a vor- 
tex tangle obeys the Kolmogorov law, the tangle is self- 
similar in the inertial range, and most energy is in the 
largest vortex. The line length density then decays as 
L ~ i~ 3 / 2 , as described in Section 4.1. If a vortex tangle 
is random and has little correlation, however, the only 
length scale is £ = L -1 / 2 , yielding a decay of L <~ t -1 , 
as described herein. Therefore, the observation of how L 
decays is helpful in understanding the structure of vortex 
tangles in QT. 

However, these studies do not explain what happens to 
QT at zero temperature at small scales. There are several 
possibilities. The first is acoustic emission at vortex re- 
conncctions. In classical fluid dynamics it is known that 
vortex reconnections cause acoustic emission. In quan- 
tum fluids, numerical simulations of the GP model show 
that acoustic emission occurs at every reconnection. 19, 2 °) 
The dissipated energy during each reconnection is ap- 
proximately 3£ times the vortex line energy per unit 
length. Wc can estimate how L decays by this mecha- 
nism. The number of reconnection events per volume per 
time is on the order of kL 5 / 2 . 16 ) Since each reconnection 
should reduce the vortex length by the order of £, the 
decay of L can be described as dL/dt = — xi£kL 5 / 2 with 
a constant xi 01 order unity. The solution is 



1 1 3 

7^72 - T3J2 + 2 Xl ^ Kt 



(24) 



In superfluid helium, especially 4 He, £ is so small that 
decay due to this mechanism should be negligible. 

The second possible mechanism is the radiation of 
sound (phonons) by the oscillatory motion of a vortex 
core. The characteristic frequencies of vortex motion on 
a scale £ are of order n/£ 2 , which is too small to cause 
effective radiation. 29 ) In order to make the radiation ef- 



fective, the vortices should form small-scale structures, 
which would be realized by two consecutive processes. 
The first process is due to vortex reconnections. 5 °) In 
a dense isotropic tangle, vortex reconnections occur re- 
peatedly at a rate of order k/£ 5 per unit volume. When 
two vortices approach one another, they twist so that 
they become locally antiparallcl, and reconnect. 8 -* They 
then separate leaving a local small-scale structure. How- 
ever, even such local structures created by single recon- 
nection events are insufficient to cause effective acoustic 
radiation. The next process occurs due to a Kelvin-wave 
cascade described in Section 5.2. The Kelvin-wave cas- 
cade could transfer the energy into the small scales in 
which the radiation of sound becomes effective. 

In actual experiments, we may have to consider the 
effects of vortex diffusion, even though these effects are 
not due to dissipation. Using the vortex filament model, 
Tsubota et al. numerically studied how an inhomoge- 
neous vortex tangle diffuses. 60 - ) The obtained diffusion of 
the line length density L(x, t) was well described by 

^jM = - X £l(x, t f + DV 2 L(x, t). (25) 

Without the second term on the right-hand side, this is 
Vinen's equation (22). 6 ) The second term represents the 
diffusion of a vortex tangle. The numerical results show 
that the diffusion constant D is approximately 0.1k. The 
dimensional argument indicates that D must be on the 
order of ft. 60 ' 

6. New QT experiments 

This section reviews briefly recent important experi- 
ments on QT, although additional research is needed for 
all of these experiments. 

6.1 Temperature- dependent transition to QT 

Krusius et al. conducted a series of experimental stud- 
ies on QT in superfluid 3 Hc -B under rotation 61-66 ) and 
observed the strong temperature-dependent transition of 
a few seed vortices to turbulence. A seed vortex can de- 
velop into turbulence when the temperature is lower than 
the onset temperature T on , whereas the seed vortex does 
not lead to turbulence above T on . The onset temperature 
T on is approximately 0.6T C with the superfluid transition 
temperature T c being independent of flow velocity. 64 ^ 
The key characteristics of using 3 He -B for research on 
QT are as follows. First, the normal fluid component is 
too viscous to become turbulent, which is very different 
from the case of superfluid 4 He. Secondly, monitoring the 
NMR absorption spectra enables us to count the number 
of quantized vortices, which is not possible in 4 He. 

Figure 4 shows schematically the vortex instability and 
turbulence in rotating 3Hc-B in the turbulent tempera- 
ture region studied by the Helsinki group. A seed vortex 
is injected into 3 He-B in a cylindrical vessel rotating with 
the angular velocity f2. In the laboratory frame, the nor- 
mal fluid component causes solid body rotation, while the 
superfluid component remains at rest. When the tem- 
perature is higher than T on , the vortex just orientates 
along the rotational axis because of the mutual friction. 
If the temperature is lower than T on , however, the vortex 
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becomes unstable, through several vortex reconnections 
that occur chiefly near the boundary, 65 ^ splitting into 
lots of vortices, and eventually reaching an equilibrium 
vortex lattice state. While the turbulent front propagates 
to the vortex-free region, the front takes on the beauti- 
ful "twisted vortex state". 66 ) The onset temperature is 
determined by considering the dynamic mutual friction 
parameter £ = (1 — a') /a, where a and a' are the mutual 
friction coefficients appearing in Eq. (6). 61 ) This param- 
eter £ works in this system like the Reynolds number in 
a usual fluid. 



n 



n 



seed vortex single-vortex turbulent 
in rotating flow instability burst 




vortex front relaxing 

and twist 
twisted state 



Fig. 4. Schematic vortex instability and turbulence in rotat- 
ing 3 He-B in the turbulent temperature region studied by the 
Helsinki group. 62 ' 



6.2 Dissipation at very low temperatures 

Recently a few experimental studies on these topics 
have been conducted, showing the reduction of the dis- 
sipation at low temperatures. Eltsov et al. studied the 
vortex front propagation into a region of vortex- free flow 
in rotating superfluid 3 Hc-B by NMR measurements fol- 
lowing a series of investigations. 67 ** The observed front 
velocity as a function of temperature shows the tran- 
sition from laminar through quasiclassical turbulent to 
quantum turbulent. The front velocity is related to the 
effective dissipation, which exhibits a peculiar reduction 
at very low temperatures below approximately 0.25T C 
with the critical temperature T c . Eltsov et al. claim that 
this is attributable to the bottleneck effect. Walmsley 
et al. made another observation of the effective viscos- 
ity of turbulence in superfluid 4 He. 30 ) Turbulence was 
produced by an impulsive spin down from an angular ve- 
locity to rest for a cube-shaped container, and the line 
length density was measured by scattering negative ions. 
The observed effective kinematic viscosity showed the 
striking reduction at low temperatures below approxi- 
mately 0.8 K. In this case, the bottleneck effect may not 
be so significant. The authors believe this may be due to 
another mechanism, namely, the difficulty in transferring 
energy through wavenumbers from the three-dimensional 
Richardson cascade to the one-dimensional Kelvin-wave 
cascade. 

6.3 QT created by vibrating structures 

Recently, vibrating structures, such as discs, spheres, 
grids, and wires, have been widely used for research 
into QT. 36 ) Despite detailed differences between the used 



structures, the experiments show some surprisingly com- 
mon phenomena. This trend started with the pioneering 
observation of QT on an oscillating microsphere by Jager 
et al. 68 ' 

The sphere used by Jager et al. had a radius of approx- 
imately 100 /im. The sphere was made from a strongly 
ferromagnetic material and had a very rough surface. 
The sphere was magnetically levitated in superfluid 4 Hc, 
and its response with respect to the alternating drive 
was observed. At low drives, the velocity response v was 
proportional to the drive Fp , taking the " laminar" form 
Fd = X(T)v, with the temperature-dependent coefficient 
A(T). At high drives, the response changed to the "tur- 
bulent" form Fd — l{T)(v 2 — Vq) above the critical ve- 
locity vq. At relatively low temperatures the transition 
from laminar to turbulent response was accompanied by 
significant hysteresis. 

Subsequently, several groups have experimentally in- 
vestigated the transition to turbulence in superfluids 4 He 
and 3 He-B by using grids, 31, 69 ~ 72 ) wires, 73 ~ 77 ) and tun- 
ing forks. 78 - ) The details of the observations are described 
in a review article. 3e ' Although a detailed discussion of 
the cases of 4 He and 3 He at both low and relatively high 
temperatures is required, we shall next briefly describe a 
few important points. 

These experimental studies reported some common 
behaviors independent of the details of the structures, 
such as the type, shape, and surface roughness. The ob- 
served critical velocities are in the range from 1 mm/s 
to approximately 200 mm/s. Since the velocity is usu- 
ally much lower than the Landau critical velocity of ap- 
proximately 50 m/s, the transition to turbulence should 
come from not intrinsic nucleation of vortices but exten- 
sion or amplification of remnant vortices. Such behavior 
is shown in the numerical simulation by the vortex fila- 
ment model. 79 ' Figure 5 shows how the remnant vortices 
that are initially attached to a sphere develop into tur- 
bulence under an oscillating flow. Such behavior must 
be related to the essence of the observations, many unre- 
solved problems remain, such as the nature of the critical 
velocity and the origin of the hysteresis for the transition 
between laminar and turbulent response. 

6.4 Visualization of QT 

There has been little direct experimental information 
about the flow in superfluid 4 He. This is mainly because 
usual flow visualization techniques are not applicable to 
superfluid. However, the situation is rapidly changing. 37 - 1 
For QT, we should seed the fluid with tracer particles in 
order to visualize the flow field and, if possible, quan- 
tized vortices, which are observable by appropriate opti- 
cal techniques. 

A significant contribution was made by Zhang and Van 
Server. 55 ' Using the particle image velocimetry (PIV) 
technique with 1.7-/im-diameter polymer particles, they 
visualized a large-scale turbulent flow both in front of 
and behind a cylinder in a counterflow in superfluid He 
at finite temperatures. In classical fluids, such large-scale 
turbulent structures are seen downstream of objects such 
as cylinders, with the structures periodically detaching to 
form a vortex street. 80 - ) In the present case of 4 He coun- 
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Fig. 5. Evolution of the vortex line near a sphere of radius 100 fim 
in an oscillating superflow of 150 mm - 1 at 200 Hz [Hannincn, 
Tsubota and Vinen: Phys. Rev. B 75 (2007) 064502, reproduced 
with permission. Copyright 2007 by the American Physical So- 
ciety]. 

tcrflow, the locations of the large-scale turbulent struc- 
tures was relatively stable, and they did not detach and 
move downstream, although local fluctuations in the tur- 
bulence were evident. 

Here, it is necessary to understand whether such tracer 
particles follow the normal flow or the superflow, or even 
a more complex flow. Poole et al. studied these problems 
theoretically and numerically to show that the situation 
changes depending on the size and mass of the tracer 
particles. 81 ^ 

Another significant contribution was the visualization 
of quantized vortices by Bewley et al. 56 -* In their ex- 
periments, the liquid helium was seeded with solid hy- 
drogen particles smaller than 2.7 /mi at a temperature 
slightly above T\ after which the fluid was cooled be- 
low T\. When the temperature is above T\, the particles 
were seen to form a homogeneous cloud that disperses 
throughout the fluid. However, on passing through T\, 
the particles coalesced into web-like structures. Bewley 
et al. suggested that these structures represent decorated 



quantized vortex lines. They reported that the vortex 
lines appear to form connections rather than remaining 
separated, and were homogeneously distributed through- 
out the fluid. The fork-like structures may indicate that 
several vortices are attached to the same particle, as indi- 
cated by numerical simulations of vortex pinning. 82 ' This 
experimental work on the visualization of quantized vor- 
tices is described in detail elsewhere. 83 - 1 

7. Summary and discussions: QT in atomic 
Bose-Einstein condensates 

In this article, we have reviewed recent research on 
QT in superfluid helium. Research on QT is currently 
one of the most important branches in low-temperature 
physics, attracting the attention of many scientists. QT 
is comprised of quantized vortices as definite elements, 
which differs greatly from conventional turbulence. Thus, 
investigation of QT may lead to a breakthrough in un- 
derstanding one of the great mysteries of nature. This 
last section addresses the topic of QT in atomic Bose- 
Einstein condensates (BECs). 

The realization of Bose-Einstein condensation in 
trapped atomic gases in 1995 has stimulated intense ex- 
perimental and theoretical activity in modern physics. 21 ) 
As proof of the existence of superfluidity, quantized vor- 
tices were created and observed in atomic BECs, and 
numerous efforts have been devoted to solving a num- 
ber of fascinating problems. 84 ) Atomic BECs have several 
advantages over superfluid helium. The most important 
point is that modern optical techniques enable us to di- 
rectly visualize quantized vortices. The long history of 
research into superfluid helium reveals two main cooper- 
ative phenomena of quantized vortices. One is a vortex 
lattice in a rotating superfluid, and the other is a vor- 
tex tangle in QT. To date, most studies of quantized 
vortices in atomic BECs have been devoted to the for- 
mer, and few have examined the latter. Recently, it has 
been shown theoretically that QT can also be created in 
trapped BECs and that the energy spectrum obeys the 
Kolmogorov law. 85 ' 86 ) Considering the advantages of this 
system, it would be possible to perform new studies on 
QT beyond those in superfluid helium. 
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